Node / meshtastic / Meshtastic-Android / files / core / model / src / commonTest / kotlin / org / meshtastic / core / model / util / AirQualityIndexTest.kt

Displaying Raw • Download

T8b949e/*

* Copyright (c) 2026 Meshtastic LLC

*

* This program is free software: you can redistribute it and/or modify

* it under the terms of the GNU General Public License as published by

* the Free Software Foundation, either version 3 of the License, or

* (at your option) any later version.

*

* This program is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the

* GNU General Public License for more details.

*

* You should have received a copy of the GNU General Public License

* along with this program. If not, see <https://www.gnu.org/licenses/>.

*/

Tff7b72package T7ee787org.meshtastic.core.model.util

Tff7b72import T7ee787kotlin.math.abs

Tff7b72import T7ee787kotlin.test.Test

Tff7b72import T7ee787kotlin.test.assertEquals

Tff7b72import T7ee787kotlin.test.assertNull

Tff7b72private Tff7b72const Tff7b72val Te6edf3SECONDS_PER_HOUR Tff7b72= T79c0ff3T79c0ff6T79c0ff0T79c0ff0L

Tff7b72private Tff7b72const Tff7b72val Te6edf3NOW Tff7b72= T79c0ff1Te6edf3_000_000L

Tff7b72private Tff7b72const Tff7b72val Te6edf3EPSILON Tff7b72= T79c0ff0.001

Tff7b72class T56d364AirQualityIndexTest Tb4b4b4{

T8b949e// Current (2024) EPA PM2.5 24-hour breakpoint table reference values (concentration -> AQI).

Tf0883e@Test

Tff7b72fun Td2a8ffpm25ToAqi_matches_epa_breakpoint_tableTb4b4b4(Tb4b4b4) Tb4b4b4{

Te6edf3assertEqualsTb4b4b4(T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff0.0Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff5T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff9.0Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff5T79c0ff1Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff9.1Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff1T79c0ff0T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff3T79c0ff5.4Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff1T79c0ff0T79c0ff1Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff3T79c0ff5.5Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff1T79c0ff5T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff5T79c0ff5.4Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff1T79c0ff5T79c0ff1Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff5T79c0ff5.5Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff0T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff1T79c0ff2T79c0ff5.4Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff0T79c0ff1Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff1T79c0ff2T79c0ff5.5Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff3T79c0ff0T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff2T79c0ff2T79c0ff5.4Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff3T79c0ff0T79c0ff1Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff2T79c0ff2T79c0ff5.5Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff5T79c0ff0T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff3T79c0ff2T79c0ff5.4Tb4b4b4)Tb4b4b4)

Tb4b4b4}

T8b949e// Regression for meshtastic/design#54: PM2.5 = 10 µg/m³ was "42, Good" on the legacy table;

T8b949e// on the 2024 table it must be "53, Moderate", matching iOS (meshtastic/Meshtastic-Apple#2075).

Tf0883e@Test

Tff7b72fun Td2a8ffpm25ToAqi_uses_2024_breakpoints_for_low_concentrationsTb4b4b4(Tb4b4b4) Tb4b4b4{

Te6edf3assertEqualsTb4b4b4(T79c0ff5T79c0ff3Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff1T79c0ff0.0Tb4b4b4)Tb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffpm25ToAqi_clamps_negative_and_above_scale_concentrationsTb4b4b4(Tb4b4b4) Tb4b4b4{

Te6edf3assertEqualsTb4b4b4(T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(Tff7b72-T79c0ff5.0Tb4b4b4)Tb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff5T79c0ff0T79c0ff0Tb4b4b4, Te6edf3AirQualityIndexTb4b4b4.Te6edf3pm25ToAqiTb4b4b4(T79c0ff1T79c0ff0T79c0ff0T79c0ff0.0Tb4b4b4)Tb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_returns_null_when_most_recent_hour_missingTb4b4b4(Tb4b4b4) Tb4b4b4{

T8b949e// Only hour 1 (an hour ago) has data - EPA requires the most recent hour (c1) to be present.

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4)

Te6edf3assertNullTb4b4b4(Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)Tb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_returns_null_with_fewer_than_two_valid_hoursTb4b4b4(Tb4b4b4) Tb4b4b4{

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4)

Te6edf3assertNullTb4b4b4(Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)Tb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_returns_null_when_second_valid_hour_is_outside_recent_threeTb4b4b4(Tb4b4b4) Tb4b4b4{

T8b949e// Two valid hours in the 12h window (now + 11h ago), but only one within the most recent 3 hours.

T8b949e// EPA requires 2 of the 3 most recent hours, so this must not report a value.

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- T79c0ff1T79c0ff1 Tff7b72* Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4)

Te6edf3assertNullTb4b4b4(Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)Tb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_averages_stable_readings_with_weight_factor_oneTb4b4b4(Tb4b4b4) Tb4b4b4{

T8b949e// No variation across hours -> weight factor stays at 1, so NowCast is a plain average.

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- T79c0ff2 Tff7b72* Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4)

Tff7b72val Te6edf3result Tff7b72= Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3resultTff7b72!!Tb4b4b4, Te6edf3EPSILONTb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_weights_recent_hours_more_heavily_when_decliningTb4b4b4(Tb4b4b4) Tb4b4b4{

T8b949e// c1=20 (now), c2=10 (1h ago). weightFactor = 1 - (20-10)/20 = 0.5 (also the EPA floor).

T8b949e// NowCast = (20*1 + 10*0.5) / (1 + 0.5) = 25 / 1.5

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff1T79c0ff0.0Tb4b4b4)

Tff7b72val Te6edf3result Tff7b72= Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff5.0 Tff7b72/ T79c0ff1.5Tb4b4b4, Te6edf3resultTff7b72!!Tb4b4b4, Te6edf3EPSILONTb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_applies_minimum_weight_factor_floorTb4b4b4(Tb4b4b4) Tb4b4b4{

T8b949e// Range far exceeds max, so the raw weight factor would go deeply negative - it must floor at 0.5.

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff1T79c0ff0T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff1.0Tb4b4b4)

Tff7b72val Te6edf3result Tff7b72= Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)

T8b949e// NowCast = (100*1 + 1*0.5) / 1.5

Te6edf3assertEqualsTb4b4b4(T79c0ff1T79c0ff0T79c0ff0.5 Tff7b72/ T79c0ff1.5Tb4b4b4, Te6edf3resultTff7b72!!Tb4b4b4, Te6edf3EPSILONTb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_averages_multiple_readings_within_the_same_hourTb4b4b4(Tb4b4b4) Tb4b4b4{

Tff7b72val Te6edf3readings Tff7b72=

Te6edf3listOfTb4b4b4(

Te6edf3NOW Te6edf3to T79c0ff1T79c0ff0.0Tb4b4b4,

Te6edf3NOW Tff7b72- T79c0ff6T79c0ff0 Te6edf3to T79c0ff3T79c0ff0.0Tb4b4b4, T8b949e// same hour bucket as above -> averages to 20.0

Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4,

Tb4b4b4)

Tff7b72val Te6edf3result Tff7b72= Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)

T8b949e// Both hourly buckets average to 20.0 -> stable, so NowCast is 20.0 regardless of weighting.

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3resultTff7b72!!Tb4b4b4, Te6edf3EPSILONTb4b4b4)

Tb4b4b4}

Tf0883e@Test

Tff7b72fun Td2a8ffcomputeNowCastPm25_ignores_readings_older_than_the_twelve_hour_windowTb4b4b4(Tb4b4b4) Tb4b4b4{

Tff7b72val Te6edf3readings Tff7b72= Te6edf3listOfTb4b4b4(Te6edf3NOW Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3NOW Tff7b72- T79c0ff1T79c0ff3 Tff7b72* Te6edf3SECONDS_PER_HOUR Te6edf3to T79c0ff5T79c0ff0T79c0ff0.0Tb4b4b4)

Tff7b72val Te6edf3result Tff7b72= Te6edf3AirQualityIndexTb4b4b4.Te6edf3computeNowCastPm25Tb4b4b4(Te6edf3readingsTb4b4b4, Te6edf3NOWTb4b4b4)

Te6edf3assertEqualsTb4b4b4(T79c0ff2T79c0ff0.0Tb4b4b4, Te6edf3resultTff7b72!!Tb4b4b4, Te6edf3EPSILONTb4b4b4)

Tb4b4b4}

Tff7b72private Tff7b72fun Td2a8ffassertEqualsTb4b4b4(Te6edf3expectedTb4b4b4: Tffa657DoubleTb4b4b4, Te6edf3actualTb4b4b4: Tffa657DoubleTb4b4b4, Te6edf3epsilonTb4b4b4: Tffa657DoubleTb4b4b4) Tb4b4b4{

Te6edf3kotlinTb4b4b4.Te6edf3testTb4b4b4.Te6edf3assertTrueTb4b4b4(Te6edf3absTb4b4b4(Te6edf3expected Tff7b72- Te6edf3actualTb4b4b4) Tff7b72< Te6edf3epsilonTb4b4b4, Ta5d6ff"Ta5d6ffexpected Tffd700$Te6edf3expectedTa5d6ff but was Tffd700$Te6edf3actualTa5d6ff"Tb4b4b4)

Tb4b4b4}

Tb4b4b4}